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1.  Introduction 


The  logistic  distribution  given  by  the  distribution  function 

F(x)  » (1  «■  exp(-(o+8x))_1  , --  < x < • , (1.1) 

where  -•<<»<•  and  0 > 0,  is  remarkably  similar  to  the  normal  distri- 
bution with  mean  » a and  standard  deviation  ■ 6ir//3.  Moreover  for  a ■ 0 
and  B ■ 1,  the  logistic  r.v.  X and  its  p.d.f.  f(-)  are  related  to 
the  d.f.  F(-)  by  analytically  simple  forms  x - log(F/(l-F))  and 
f(x)  » F(x)(l-F(x)) . Because  of  these  reasons  this  distribution  is  widely 
employed  as  a substitute  for  the  normal  distribution  in  applications  such 
as  bioassay  and  quantal  response  data  problems,  e.g.  Berkson  (1944),  and 
growth  curve  analysis  (Verhulst  (1845)).  The  importance  of  the  logistic 
law  in  modelling  of  stochastic  phenomena  has  resulted  in  numerous  studies 
involving  probabilistic  and  statistical  aspects  of  this  distribution.  For 
example,  Gumbel  (1944),  Gumbel  and  Keeney  (1950),  and  Talacko  (1956)  show 
that  it  arises  as  a limiting  distribution  in  various  situations;  Birabaum 
and  Dudman  (1963),  Gupta  and  Shah  (1967)  study  order  statistics  from  it; 
and  many,  e.g.  Antle  et  al.  ( 1970)  and  Tartar  and  Clark  (1965)  study  in- 
ference questions  concerning  it.  For  an  excellent  summary  of  results  about 
logistic  distribution  refer  to  Chapter  22  of  Johnson  and  Kotz  (1970). 

As  might  be  expected  in  view  of  the  similarity  between  the  logistic 
and  normal  distributions,  the  sample  mean  and  the  sample  variance,  the 
moment  estimators  of  a and  8,  are  effective  tools  for  statistical  deci- 
sions involving  the  logistic  distribution.  Antle  et  al.  (1970)  give  a func- 
tion of  the  mean  as  a confidence  interval  estimate  of  a when  B is  known. 


Schaffer  and  Sheffield  (1973)  make  a case  that  in  terns  of  mean  squared 
error  the  moment  estimators  are  as  good  estimators  of  a,  6 as  their  M.L.E.'s. 

The  fact  that  the  distribution  of  the  sample  mean  has  monotone  likelihood  ratio  with 
respect  to  a when  0 is  known  leads  to  a uniformly  most  accurate  confi- 
dence interval  for  a and  UMP  tests  for  one-sided  hypotheses  involving  a. 

The  sampling  distribution  of  the  mean  is  a primary  requirement  for  these 
statistical  purposes.  The  papers  due  to  Antle  et  al.  (1970)  and  Schafer 
and  Clark  (1965)  both  contain  Monte  Carlo* ed  results  for  this  distribution. 

Goel  (1975)  obtains  an  expression  for  the  d.f.  of  the  sum  of  i.i.d.  logistic 
variates  by  using  the  Laplace  transform  inversion  method  for  convolutions 
of  Polya  type  function  developed  by  Schoenberg  (1953)  and  Hirschman  and 
Widder  (1955).  He  uses  this  approach  because  an  inversion  of  the  character- 
istic function  of  the  sum  leads  him  to  a very  slowly  converging  series  for 
its  d.f.  In  this  paper  we  demonstrate  that  the  characteristic  function  can 
be  inverted  to  obtain  a neat  expression  for  the  d.f.  Also  by  noting  the 
length  of  the  tails  of  the  distribution  we  develop  a simple  student-t  approxi- 
mation for  it  and  compare  it  with  the  normal  and  an  Edgeworth  series  approxi- 


2.  D.F.  of  the  Sum  of  I.I.D.  Logistic  Variates 


2 - I (2.2) 

i-l  1 

Goel  (1975)  obtains  the  expression 

Fn(x)  - 0.5  ♦ i l ((m+n-l) l/m!)Re((c*id)n)  , (2.3) 

where  c • (■  ♦ n/2)/a,  d - x(s/3)**2a,  a - (m  ♦ n/2)2  ♦ (x2  n/12)  for  the 
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d.f.  of  Z by  applying  the  inversion  formula  directly.  As  he  notes,  this 

series  converges  slowly  and  is  therefore  computationally  inconvenient.  In 

this  section  we  use  the  inversion  formula  for  characteristic  functions  to 

obtain  a finite  series  expression  for  ?„(*)•  The  method  used  requires 

only  a knowledge  of  contour  integration  and  the  Mittag-Leffler  expansion 

of  the  characteristic  function  of  Z which  may  be  found  in  sources  such  as 

the  monograph  by  Widder  (1972). 

The  characteristic  function  +n(*)  of  Z is  given  by 

4 (s)  » (iris/sin  iris)n.  (2.4) 

n 

Its  Mittag-Leffler  expansions  are  given  by 

(iris/sin  iris)"  ■ (is)n  £ [ A . /(is+k)n'p,  (2.5) 

k»—  p«o  n_1,p 

when  n is  even,  and 

(iris/sin  iris)n  - (is)n  l l (-l)kA  , /(is*k)n"p  (2.6) 

k-—  p-o  n"1,P 

when  n is  odd. 


The  constants  AR  in  (2.5)  and  (2.6)  can  be  computed  by  using  the  follow- 
ing equations:  / 


(irs/sin  irs)n  ■ A„  ♦ A_  . s ♦ ...  ♦ A„  sn, 
n9u  n9i  n9n 


(2.7) 


and 


ws/sin  irs  • l (-l)k_1(22k-2)/(2k)l  B,.  (its)". 


1 2k 


k-0 


2k' 


(2.8) 


the  B2k's  being  Bernoulli  numbers . Since  the  d.f.  of  Z is  symmetric 
t zt 

gives,  0 < x < b < • , 


about  zero,  it  suffices  to  find  F (x)  for  x > 0.  The  inversion  formula 

n 


, N+%  -isx  -isb  n 

».(«•».«  - TT  < b <*•« 


Consider  the  contour  integral 

■ixz  _-ibz> 

I,  ‘ '' 


(2.10) 


. ji/  a . 

Si 

where  is  a clockwisely  oriented  contour  given  by  CN  ■ ♦ CJJ  . 

» {z  : | z | < and  Im  z ■ 0} 

and  (2.11) 

CJj  » {z  : |z|  ■ N^*s  »nd  Im  z < 0)  . 
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It  is  clear  that 


provided  that 


F_(b)-F_(x)  ■ lim  IN(b,x) , 

n n N-m*  n 


i , u-lxl  - .-ibzi,  .t« 
v, 27  /,  ^ li  ^sin  iriz 


■)  dz  » 


(2.12) 


(2.13) 


(2.13)  is  easy  to  establish  because 

If  t'ixz  , It  ,n.T  . , 2i(N.lj)  .rle 

\i  —n~  sin  it> dl  i c7nR5T77TTO)  LI  - 


■ixz  | 

■jj—  dz  ♦ 0 (2.14) 


N 

as  N Now  we  use  (2.4)  to  evaluate  F (x)  when  n is  even: 

i n_1  • / -ixz  .-ibzx  n 1 

Ub.x)  - U/  l l AnDi^ ^_L(i2)n-l 

N 2 CN  p-0  k—  n,p  (iz*k)n_p 

n-1  N _ , _ . / -ixz  -ibzt 

• l l \ n iP  57T  / 2 =4,  1 <*2 

p-0  k-1  n,p  (Zfik)n'p 

n-l  . " n n 1 (»~ixz  .-ibz\ 

. T ‘ / ; A i*  t"'1  il ^ L dz 

p.o  171  ‘ k-N.l  "•»  (t.ik)"-P 

n-l  . Si  - _ _ - / -ixz  -ibzx 

♦ l 1 / l a ipzn_1^ ^2 Uz. 

pio^Lkio  n*p  (z-ik)n-p 


(2.15) 


(2.16) 


The  integrand  of  equation  (2. IS)  has  a removable  singularity  at 
and  poles  of  order  (n  -p)  at  z ■ -ik,  k - 1,2,...  . Hence  the  inti 

grands  in  the  last  two  terms  of  equation  (2.16)  are  analytic  in  C^.  There- 




zero. 


Now  let 


n-l  N , . / -ixz  -ibzv 

Ub.X  - £ I *nniP5Ir/  1 15 ■>  dt. 

N p-Okil  "-P  (l.ik)"'P 

N - n-l  -iyz 

IN(y)  ■ l -r^r  J jpr-dz. 

N k-1  2iri  (z*ik)n’P 

N 

I«(y)  ■ I Res(-ik), 

" k-1 


(2.17) 


(2.18) 


n-l  -iyz 

where  Res(-ik)  is  the  residue  of  5 =— — - at  z ■ -ik,  k - 1,...,N 

(z*ik)n_p 

and  is  given  by 


• (-iiph  gi!p  x-1> 


(2.19) 


QZ  ! X— ik 

- Turwaiire^  c. 

r-0  p r r 
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Thus 


n-l-p 


r » 


1U  I„(y)  - I (-l)r*1(-i)P  <";})  C l kP*r  .-‘y.  (2.21) 

r«C  v ' k«l 


In  the  development  from  (2.12)  if  we  let  b -*■  • then 

1 - F (x)  - lim  lim  IN(b,x) 

n , _ n 


X X*  (-1,r41  Vp  <£>  £ J.  kP*r 


-kx 


p*0 


Using  induction,  the  infinite  series  \ k**+r  e”kx  can  be  shown  to  con- 
p-*r  k«l 

verge  to  £ s (m-l) ! (e'*/(l-e  ))  , where  s 's  are  Stirling's 

sal  P**»® 

number  of  the  second  kind  given  by 

s».r.  * Ctp-r)!)'1  I (-l)1  (*K»-i)P*r.  (2.22) 

i»0 

By  a similar  method,  an  expression  can  be  obtained  for  Fn(*)  when 
n is  odd.  These  results  are  stated  in  the  following  theorem* 


THEOREM  2.1.  Let  Xj,  . . . , Xfl  be  i.i.d.  with  d.f.  given  by 

n 

(2.1).  Then  the  d.f.  of  Z«  [ X.  is  given  by 

• * 


i-1 

n-1  n-p-1  p-*>r*l 


-kx 


1 • F»“>  • Jo  Jo  J, 

(2.23) 

kx 


when  n is  even,  and  by 

n;l  iwp-1  p.r*l  (_i<r.k.i 
p»0  r«0  k«l 


n-l  n-p-l  p+ 

1 - Fn(z)  - I l l (-D™*  A_  _s. 


»,p>r.l.k  JT  - 


when  n is  odd. 


(1*0 

(2.24) 


xvlT  * 


3.  Some  /^proximations  for  the  Convolution  of  Logistic  R.V.'s 

As  mentioned  earlier  the  logistic  distribution  is  very  similar  to  the 
normal  distribution.  It  is  therefore  to  be  expected  that  the  normal  distri 


Ik 

■ *. 
h* 
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bution  would  reasonably  approximate  the  distribution  of  the  sum  of  k in- 
dependent logistic  r.v.'s  for  k >_  2.  Specifically,  if 

1 - ^3 7F  Z/tt  (3.1) 

denotes  the  standardized  sum  of  k i.i.d.  logistic  r.v.'s;  then  we  have  a 

simple  approximation: 

* *U1.  (3.2) 

* * 

where  Fk(*)  is  the  d.f.  of  Z and  « is  the  standard  normal  d.f. 

This  approximation  can  be  improved  by  applying  a two  term  Edgeworth  correc- 
tion to  it  to  obtain 

F*(z)  * *(z)  - C20k)'1(z3-3z)w(.z),  (3.3) 

where  W(z)  ■ ( 2tt ) -Js  exp(-z2/2). 

Although  the  above  approximations  are  good,  it  may  be  possible  to  im- 
prove on  them  by  using  student's  t distributions  which  is  similar  to  the 
normal  distribution  in  shape  but  has  relatively  long  tails.  The  degree  of 
freedom  v of  the  approximating  t distribution  can  be  obtained  by  equating 


the  coefficients  of  kurtosis. 

s2(z)  - 

* 

(3.4) 

which  gives 

v * 

Sk  ♦ 4 

• 

(3.5) 

Thus  considering  the 

standardized  Z 

and  t. 

★ * 
i.e.  Z and  ty. 

respec- 

tively,  where 

• 

t ■ 

V 

t /a(t  ) 

V ' V 

« /(v 

-2 }/\>  t 

V J 

(3.6) 

this  approximation  is 

given 

by 

* P“5k.4i2/ 

/ SkTT  . 

5FT  K 

(3.7) 
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4. 

Figure  1 and  Table  1 illustrate  the  quality  of  the  three  viz.  the 
normal,  the  Edgeworth-corrected  normal,  and  the  student's  t,  approximations 
for  k = 2 and  3,  respectively.  The  maximum  error  for  the  t -approximation 
for  k * 2 is  about  0.005  and  decreases  to  0.0007  when  k * 3.  Goel  ( 1975 ) 
gives  a table  of  cumulative  d.f.  of  the  standarized  mean  of  samples  from 
logistic  population  for  k ■ 10  correct  to  seven  decimal  places.  We  believe 
that  for  this  value  of  k the  t -approximation  would  yield  results  of  compar- 
able accuracy. 


Evaluation  of  the 


jroximation 
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TABLE  1. 


A Comparison  of  Three  Approximations 
k - 3. 


X 

f;u) 

F*(x)-*(x) 

F*(x)-G2(x) 

F*(x)-H3(x) 

0.05 

0.5209 

0.0010 

0.0000 

0.0001 

0.15 

0.5625 

0.0029 

0.0000 

0.0003 

0.25 

0.6033 

0.0046 

0.0008 

0.0005 

0.45 

0.6809 

0.0073 

-0.0017 

0.0007 

0.65 

0.7506 

0.0084 

-0.0006 

0.0007 

0.85 

0.8106 

0.0083 

-0.0007 

0.0007 

1.00 

0.8486 

0.0073 

-0.0008 

0.0004 

1.20 

0.8903 

0.0054 

-0.0007 

0.0002 

1.45 

0.9291 

0.0026 

-0.0004 

0.0000 

1.75 

0.9598 

-0.0001 

0.0001 

-0.0002 

2.50 

0.9918 

-0.0020 

0.0004 

-0.0002 

3.00 

0.9975 

-0.0012 

0.0001 

0.0000 

F.(x)  * Distribution  function  of  the  standardized  sum  of  3 i.i.d. 

logistic  r.v. 's 

♦ (x)  ■ Distribution  function  of  the  standard  normal  r.v. 

G^(x)  » Edgeworth  series  approximation  given  by  right  hand  side  of 

equation  (3.3). 

H3(x)  - 


Distribution  function  of  standardized  student's  t r.v.  with 
19  degrees  of  freedom,  as  given  by  equation  (3.7). 
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